Non-equilibrium Green's function approach to inhomogeneous quantum many-body 
systems using the Generalized Kadanoff Baym Ansatz 
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In non-equilibrium Green's function calculations the use of the Generalized Kadanoff-Baym 
Ansatz (GKBA) allows for a simple approximate reconstruction of the two-time Green's function 
from its time-diagonal value. With this a drastic reduction of the computational needs is achieved in 
time-dependent calculations, making longer time propagation possible and more complex systems 
accessible. This paper gives credit to the GKBA that was introduced 25 years ago. After a detailed 
derivation of the GKBA, we recall its application to homogeneous systems and show how to extend 
it to strongly correlated, inhomogeneous systems. As a proof of concept, we present results for a 2- 
electron quantum well, where the correct treatment of the correlated electron dynamics is crucial for 
the correct description of the equilibrium and dynamic properties. 

PACS numbers: 05.10.-a, 05.30.-d, 71.10.-w 



I. INTRODUCTION 

For the time-dependent description of non- 
equilibrium processes the method of non-equilibrium 
Green's function (NEGF) has been widely used, since it 
allows for self-consistent treatment of electron-electron 
correlations, non-perturbative inclusion of external 
fields and systematic approximations via Feynman 
diagrams. The central quantity is the one-particle, 
two-time function, G(t,t'), the time-evolution of 
which is governed by the Keldysh-Kadanoff-Baym- 
equations QJ. To numerically solve these equations a 
self-energy is introduced, which can be determined by 
many-body perturbation theory and which leads to a 
closed equation for G. Still, as G inherently depends on 
two times t and t' , the time propagation is numerically 
demanding and the memory and CPU time needs scale 
quadratically with the propagation length, see e.g. 121 . 

This restriction can be drastically alleviated by the 
introduction of a further approximation, the general- 
ized Kadanoff-Baym ansatz (GKBA), which was intro- 
duced by Lipavsky, Spicka and Velicky some 25 years 
ago 1 3 1. With the GKBA, for each time propagation step, 
the two-time Green's function is reconstructed from its 
time-diagonal value: G(t,t') = F GK ^ A [G(t = t')]. Asa 
great advantage it reduces the amount of needed mem- 
ory to a linear scaling with propagation length, since 
for the determination of G(T, T") for time-arguments 
T < T only the knowledge of G(t = t') for all t < T 
is sufficient. This simplification has made numerous ap- 
plications for spatially homogeneous systems possible. 
Here we demonstrate that the GKBA may be equally 
successful in computing the behaviour of finite inhomo- 
geneous systems. 

The paper is organized as follows: After a short rec- 
ollection of the basics of the NEGF formalism we, in de- 
tail, derive and list the properties of the GKBA. There- 



after we give a brief overview about its application to 
homogeneous systems in different fields of physics. In a 
third part we extend the GKBA to inhomogeneous sys- 
tems, using the technique of adiabatic switching (AS) |4J 
to correctly obtain the associated correlated initial state. 
Finally the applicability of the GKBA to the spctrum of 
a two-electron quantum well is tested at different cou- 
pling strengths. 



II. THEORY 

A. Non-equilibrium Green's Functions (NEGF) 

To describe correlation effects and excitations in quan- 
tum many-particle systems we chose the NEGF ap- 
proach, as it allows for a systematical inclusion of cor- 
relations by diagrammatic expansions. In contrast to 
density matrix based schemes, the Green's function 
method additionally easily offers direct access to dy- 
namical spectral information as well as particle removal 
and addition energies. The main quantity is the one- 
particle Green's function, defined as (we set h = 1) 



G(t,t') = s(Tc[nmHt')]), 



a) 
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where the brackets denote thermodynamical averaging 
and 7c is the time ordering operator on the Schwinger- 
Keldysh contour C |5|, on which t and t' are defined. 
\E ,( ^) denotes a one-particle annihilation (creation) oper- 
ator in a one-particle basis in second quantization. The 
equations of motion for G are the Keldysh-Kadanoff- 
Baym equations (KBE) 

{id t -h{t)]G(t,t') = 8{t-t l ) (2) 

+ J dtW{t,t)G {2) {tt-t't + ) 1 

[-ify - h{t')]G{t,t') = S(t -t') 

diW(t',t)G^{ti; t'i+), 
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where h denotes the one-particle Hamiltonian, 
G^(tt';t[t' 2 ) = -(7c [*(t)*(t / )* t (^)* t (ti)J) is the 
two-particle Green's function and W is an arbitrary in- 
teraction potential. The KBE are the first equations of 
the Martin-Schwinger Hierarchy |6|, which describes 
the coupling of the evolution of the one-particle Green's 
function to the two-particle Green's function, which it- 
self is coupled to the three-particle Green's function 
by a similar equation. To decouple this hierarchy and 
to make the KBE numerically tractable, a self-energy 
T,(t,t') = £[G(t,i')] is introduced. This self-energy 
can be found from a diagrammatic expansion in terms 
of Feynman diagrams, where only some classes of di- 
agrams are chosen according to the properties of the 
examined system. With this the KBE attain a formally 
closed form: 

[8d t -h(t)]G{t,t') = S(t-t')+ { dtE[G](t,t)G(i,t') , 



[-Qdt, - h(t')] G(t, t') = S(t - t') + J dtG(t,F)E[G](t,lf) . 

(3) 



B. Reconstruction Problem 

1. Keldysh representation and Dyson equation 

The two-time structure of the time contour C suggests 
the use of a matrix representation for G according to 
the different time orderings. There exist different repre- 
sentations, that are connected by a Keldysh rotation [7j. 
Here, we use the set involving less-, retarded and ad- 
vanced Green's functions (G < , G R and G A ) according to 
Langreth and Willkins 0, 



G 




(4) 



Note, that this representation implies a two-time depen- 
dence of each component and is overcomplete as the 
conjugation relation [G R (t,t')]* = G A (t, if) holds. To 
simplify the notation we will make use of the greater 
Green's function G > , which relates to the other compo- 
nents as: 

G>(t,t') = G R (t,t') - G A (t,t') + G<{t,t') . (5) 

With these definitions we arrive at a formal solution for 
G by time integration of Eq. <j3j yielding the nonequilib- 
rium version of the Dyson equation, 

G = G + G SG. (6) 

Here G denotes the non-interacting or Hartree-Fock 
(HF) Green's function, whose inverse is given by Gq~ : 



G 1 — G 1 • 1 — 



G n 







G n 



(7) 



where (Go) l , the common inverse of the components 
G R and G A , reads 



(Go)- 1 (t,t , ) = 5(t,t , )[W t ~h(t)- 



(8) 



The matrix multiplication in Eq. |[6} is to be understood 
as also including a time integration on the contour over 
intermediate time coordinates, so that, e.g. 

(G SG)(i',t) = J <£tZtG (t,i)S(t,!)G(i,t'). (9) 



2. Equation of motion for G < in terms of the density matrix 

Following V. Spicka et al. |9|, an intermediate step to- 
wards the GKBA is to express the equation of motion 
of G < in terms of the density matrix p(t) = — iG < (t, t). 
To start with, we provide some useful relations between 
qR/</a anc j tne reS p ec tive self-energies T, R ^ < ^ A . 

By right-multiplication of Eq. ^ with G _1 and left- 
multiplication with (Go) -1 we attain, taking the re- 
tarded/ advanced component: 



(Go) 



Q-l,R/A + Y, R / A . 



(10) 



Now taking the less-component of Eq. |6}, after left- 
multiplication with (Go) w e find 



((Go)- x G) 



(EG) 



Using the Langreth- Wilkins Rules |8| it follows: 



(Go)- 1 ^'G < + (G ) 



-i,< 



G = STG 



£<G 



(11) 



(12) 



Note, that the multiplication is to be understood in the 
same manner as in Eq. Q, including contour time inte- 
gration. Since (G ) 1:< = 0, Eq. 1 12 simplifies to: 



(G )~ 1 G < 



S fl G< +E<G 



<r< A 



and use of Eq. 1 10 1 yields: 



(G 



-1,R 



G< = EG 



Rn< 



<r<A 



E < G 



(13) 



(14) 



Analogously we find the conjugate equation, resulting 
in two final differential equations: 



G~ 1 ' R G < = E < G A , 
G < G- 1 - A = G fl S < . 



(15) 
(16) 



Now it is convenient to also split G < into two parts cor- 
responding to the time arguments t > t' and t <t'\ 

G < {t,t') = G<(t,t')-G<{t 1 t'), (17) 
G<(t,t') = Q(t-t')G < (t,t'), (18) 
Gi(t,t') = -G(t'-t)G<(t,t'). (19) 
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This allows us to separately derive an equation for G^ 
and G^ from which Eq. |l7| allows us to recover the 
equation for G < . For G^, one calculates: 



[G-^G<] (t,f) 

= [((Go)- 1 -!]^ G<] (t,f) 

dt [5(t, t) (Qdt - hit) ) Q(i- t')G < (i, t')] 



(20) 



-J di£ fl (M~)G<(M') 



dt [5(t, t) (QS(t, t')G < (t, t') + B0(t - i')%G < (*, f) 



-h(t)e(i-t')G < (i,t'))] - J dfS R (t,t)G<(f,t') 
5(t-t')fiG < (t,t') 

dt <y(t, t)e(t - o (fia t - - /*(*)) g< {t, t') 



dtY, R (t,t)G<{t,t') 



: S(t^t')QG < (t 7 t r ) 

+ B(t-t') { dtS(i - t')(G y 1 (t,t)G < (t,t') 



-J diZ R (t,t)G^(t,t') 
- 5{t-t')m < {t,t') 



+ e(t-t')( J dte(i-t , )G- 1 ' R (t,t)G < (t,t / ) 



+ [ dtZ R (t,t)G<(t,t')- [ dtT, R {t,t)G 
Jc " Jc 



Rt* ^-<(i,t r 



where the inclusion of the last term under the 0- 
function is justified, since the contour product of two 
retarded functions is again a retarded function. The two 
last terms cancel and one finds, employing Eq. p5|, 



5{t-t')m < {t,t l ) 



(21) 



+ e(t-t') I die{t-t')G- hR it,i)G < {i,t') 

S(t-t / )QG < (t,t') 

+ e(t-t') J diT, < (t,t)G A (t,t') 

-e(t'-t) / die(t'-t)G- 1 ' R (t,t)G < {t,t'). 

Jc 



Here, the last term compensates the step function in the 
second term. With Eq. ([l0| and, noting that 



dt Q(t' - t) (G ) _1 (t, t)G< (t, t') = 0. 



it follows 

S(t-t')8G < (t,t') 

+ <d{t-t r ) J dtZ < {t,t)G A (t,t / ) 

-Q(t'-t) J diQ(t' ~i)G- 1 ' R {t,i)G < {i,t r ) 
= 6(t-t')6G < (t,t') 



(22) 



+ Q(t - 1 



') J dte(t' -t)T l < {t,i)G A (i 1 t') 

+ Q(t-t') [ dte(t'-t)E il (i,£)G < (t,t')» 
Jc 

where, in the second term, the fact was used, that G A 
includes a step function by definition. Left-multiplying 
by G R and explicitly writing out the integrals, one ar- 
rives at the equation for G^(t,t') = G < (t,t / ), valid for 
times t > t' > to- 

G<{t,t') = -G R {t,t')p{t') (23) 
+ f dt I diG R {t,i)J: < {lt)G A (t,t') 

Jt' Jt 

dt / dtG R (t,t)Y, R (t,t)G < (i,t') . 

J t 

In a similar manner one can derive the equation for 

G < (t,t / ) = -G^(t,t') in the time domain t < t < t', 
which reads: 

G<(t,t') = p(t)G A (t,t') (24) 
+ I dt I dtG R (tJ)T l < (t,t)G A (t 7 t') 

Jt' Jt 

+ [ dt I diG < {t ) t)Y> A {t,t)G A {i,t') 1 

Jt 1 Jt 

Note, that by exchanging (<<^>>) and replacing the den- 
sity matrix p = : / < by/ > = 1 — / < in equations 1 23 1 and 
1 24 1, the analogous expression for G > is easily obtained. 



C. The Generalized Kadanoff-Baym Ansatz 



Combining equations d23) and 1 24 1, retaining only the 
terms without the integrals^ the Generalized Kadanoff- 
Baym Ansatz |3] is recovered: 



gIkba(M') = -G K (M')/«(0 + /«WG A (M')- (25) 

With this equation a means for the reconstruction of the 
off-diagonal Green's function from the density matrix is 
found, which obeys particle number conservation, has a 
per se causal structure and does not depend on assump- 
tions about near equilibrium, for a recent discussion, see 
eg. Ref. EJ. 
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1. Choice of propagators G K I A 

When looking at Eq. d25|, it should be noted, that 
the GKBA is only formally closed in terms of p, since 
it depends on the knowledge of the retarded (ad- 
vanced) propagators G R (G A ), which themselves obey 
non-Markovian two-time equations of motion of simi- 
lar complexity. This can be overcome by the use of prop- 
agators, which obey a Markovian evolution. In this pa- 
per we choose Hartree-Fock propagators, as they incor- 
porate the interaction at mean-field level in contrast to, 

e.g., ideal propagators. The HF propagators G^/p are 
given by: 

G^ / /(t,*') = = Ffi0[±(*-*')]expf-fi^ dffl"(*)V (26) 

where H denotes the mean-field HF Hamiltonian, 
which is governed by the time-dependent density ma- 
trix. 



2. Choice of self-energy 

To show the advantages of the GKBA we exemplarily 
apply the second order Born approximation, leading to 
a self-energy E2s(t, f), which is given by 

x< B (t,t')=s(t-t')x< F (t) 

+ G< (i, l!)W{t)W{l/)G > (f, t)G< (t, t') 

- G< (t, t')W(t)G> {t', t)G< (t, t')W{t') . 

(27) 

The evaluation of the collision integral I = J c SG, the 
right hand side of the KB equations Eq. j3j, is twofold 
simplified by the GKBA as we can use the reconstructed 
G(t, t') in / as well as in E. For instance, the less-part of 
Tj2b for t > t' now reads: 

Z< B (t,t')=5(t-t')Z< F (t) 
+ G R (t.t')f<(t') 

[W{t)W{t')f>{t')G A {t',t)G R {t,t')f<{t') 

- W{t)f>{t')G A {t\t)G R {t,t')f<{t')W{t')] . 

(28) 

So only the single-time quantity /< has to be stored in 
memory, as the HF propagators G R / A (t, t') can be com- 
puted each time step, and the demand scales linearly 
with propagation length. That way, in total three ap- 
proximations were taken: 

i. The self-energy was introduced to decouple the 
Martin-Schwinger hierarchy. 

ii. The two-time Green's function was reconstructed 
within the GKBA according to Eq. | [25| , neglecting 
the integral terms in equations |23) and <[24|. 



ill. The propagators G R I A were approximated by HF 
propagators, which can be explicitly calculated for 
each (t,t'). 



III. APPLICATION TO HOMOGENEOUS SYSTEMS - 
OVERVIEW 

The GKBA has been applied to spatially homoge- 
neous charged particle systems already in the mid 1990s 
and has allowed for many numerical studies of ultra- 
fast carrier relaxation. Electron-phonon scattering in 
semiconductors was investigated e.g. in Refs. IfTOl [Til . 
Electron-electron in plasmas and semiconductors as 
well as tests of the GKBA against full two-time calcu- 
lations were performed in Refs. H2TIT4I . The use of cor- 
related spectral functions (beyond Hartree-Fock propa- 
gators) was analyzed in Ref. lTl5l . For an overview on 
the GKBA and applications to electron-hole plasmas in 
semiconductors, see the text book 1116] . The second type 
of applications was devoted to dense plasmas, in partic- 
ular laser plasmas. To capture strong field effects in the 
Coulomb scattering process (such as harmonics genera- 
tion, inverse bremsstrahlung), a gauge-invariant gener- 
alization of the GKBA was derived in Ref. [17] and used 

in nana. 

Since application of the GKBA transforms the NEGF 
approach into a single-time theory, there should ex- 
ist close connections to the purely single-time ap- 
proach of reduced density operators (quantum BBGKY- 
hierarchy). In fact, this has been studied in detail in Ref. 
l20l , and a one to one correspondence can be established 
in the limit of free or Hartree-Fock propagators. There 
it was also shown that the GKBA does not destroy the 
conserving properties of the underlying NEGF approx- 
imation as long as no appoximations to the time struc- 
ture of the propagators are introduced. In particular, in 
the relevant case of free or HF propagators the conser- 
vation laws and sum rules are preserved |21j. For a re- 
cent overview on the GKBA and its relation to quantum 
transport and density functional theory, see [22 J. 



IV. APPLICATION TO INHOMOGENEOUS SYSTEMS 
A. Initial state preparation under the GKBA 

In order to solve the KBE Eq. ||3}, one has to supply an 
initial value for G(t ,to). For large, homogeneous sys- 
tems described in Section [ill] the initial state is reached 
from an arbitrary state after a characteristic equilibration 
time. For small, strongly correlated systems, no equili- 
bration takes place and the preparation of the correlated 
initial state has to be performed differently to avoid 
strong artifacts, such as unphysical oscillations. For full 
two-time calculations, the initial state can be found by 
means of the extended Matsubara-Schwinger-Keldysh 
time-contour (containing an additional imaginary time 
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Time t [arb.u.] 



40 45 50 



Figure 1. Fermi function for adiabatic switching. The half- 
time is to = 25 and the switching time constant is chosen to be 

r = 3 



branch) and solution of the Dyson equation Eq. |6} on its 
imaginary branch, e.g. |23|, using the same self^energy 
(for details and implementation see e.g. l24l ). Here, 
however, the application of the GKBA poses a prob- 
lem, since no corresponding equilibrium approximation 
is known so far. This can be remedied by using adiabatic 
switching as discussed below. 

In general, according to the Adiabatic Theorem 
the ground-state of an interacting system can be found 
by taking the ideal system, for which the ground-state 
is assumed to be known, and adiabatically switching on 
the interaction. If this procedure is performed slowly 
enough, the system is undergoing a transition through 
successive eigenstates of the respective Hamiltonians 
with gradually higher interaction strength. Finally, the 
fully interacting ground-state is reached. Here, we use 
this method for generating the initial state that is con- 
sistent with the GKBA. We use a monotonic switching 
function /, which has the following properties: 

/(to)=0, /(*/)= 1, 0</(i <*<*/) <1, (29) 

where tt is the end of the switching process and the start 
of the time-dependent calculation. It is important that, 
with this method, the interaction that normally is time- 
independent, becomes inherently time-dependent. This 
especially needs to be kept in mind, when dealing with 
quantities involving interaction terms. For example, the 
HF propagators now include W(t) and the self-energy 
contains Wit) and W(tf). In the calculations below we 
use a Fermi-like switching function f[ o (t) (see Fig. [TJ: 



fL(t) = i 



l+exp( t ^)' 



(30) 



where fj (to) = \ and the smoothness of the switching 
increases with the value of time constant r. 



B. Application to electrons in quantum wells 

To test the ability of the GKBA to describe correla- 
tion effects in an inhomogeneous system, we study two 
electrons in a quantum well potential. After preparing 
the system in the correlated initial state via adiabatic 
switching the system is disturbed by a short dipole kick 
EoTl with sufficiently small amplitude, and the time- 
dependent dipole moment is computed. Fourier trans- 
formation then yields the correlated dipole excitation 
spectrum in linear response with the relevant vertex cor- 
rections, thereby fully preserving conservation laws and 
sum rules |26|. All calculations are performed in the 
context of a FE-DVR basis, which drastically reduces the 
numerical complexity, for details the reader is referred 
to refs. Iffll27l 

We consider N = 2 electrons in a quantum well 
potential, which is effectively a one-dimensional sys- 
tem, if the lateral electronic motion is neglected. We 
assume, that the system is in singlet configuration 
\S,Ms) = |0, 0) . The confinement energy is given 
by Eq = h 2 /{m* L 2 ), where m* is the effective mass 
of the electrons within the quantum well potential 
of width L. The 2-particle Hamiltonian in units of 

E* reads H 2 (t) = (£?=i + f ■ Xi S(t - t )) + 

A* [(xi — x 2 ) 2 + k 2 ] , where the particle po- 

sitions. The first term denotes the single-particle con- 
tributions involving the potential energy and the dipole 
delta-kick with amplitude /q -C 1. 

The second term describes the two-particle Coulomb 
interaction. A cut-off parameter k has been inserted 
to regularize the ID-Coulomb potential, which is set 
to a value of 1 throughout the calculations. The rela- 
tive interaction strength between the electrons is given 
by the dimensionless coupling parameter A* = L/a^ = 
e 2 m*L/(4TT€oe*h 2 ), where e* denotes the material's di- 
electric constant that enters the effective Bohr radius aj$ . 



C. Numerical results 



In Fig. [2] the ground state dipole excitation spectrum 
of the 2-electron quantum well is presented for differ- 
ent values of the coupling parameter A. In gray-scale 
the results from GKBA calculations using second order 
Born self-energy are shown. The white dashed lines 
represent the excitation energies from exact diagonaliza- 
tion (ED). Exemplarily for A* = 1, the respective time- 
dependent dipole moment d(t) can be seen in the figure 
above the spectrum. The first 50 a.u. of the propagation, 
where the dipole moment is zero, accounts for the adia- 
batic switching (denoted "AS" in the figure). While d(t) 
appears to be monochromatic, in fact, it contains numer- 
ous additional frequencies which can only be resolved 
using a sufficiently long time propagation. The present 
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Figure 2. Dipole excitation of a 2-electron quantum well. 
Top: Time-dependent dipole moment for coupling A* = 1 
for GKBA with adiabatic switching ("AS"). Bottom: Dipole 
excitation spectrum [Fourier transform of d(t)] for different 
coupling strengths A*. GKBA results with second Born self- 
energy (gray-scale) are compared to the results from exact di- 
agonalization (dashed white lines). The letters (numbers) re- 
fer to one-(two-)electron excitations. Note, that the effective 
energy Eq scales oc (A*) - , so the excitations nearly constant 
with A* in this figure, effectively also scale oc (A*)~ 2 . 



le+06 




0.01 I . — . — . — 

1 10 100 1000 

time steps T 

Figure 3. Performance of second order Born NEGF simulations 
for a two-electron quantum well. GKBA lower (blue) lines vs. 
full two-time propagation, upper (red) lines. The total elapsed 
computation time is shown as a function of the number of 
propagation time steps T for fixed step size St = 0.01 a. it.. 
The dashed lines are guides to the eye indicating the expected 
asymptotic behaviour. 



GKBA calculation makes this possible. It has a total du- 
ration of T = 40000 time steps and is readily performed 
within 24 hours for a few tens of basis functions. 

Let us first discuss the general structure of the excita- 
tion spectrum that is obtained from the ED calculations. 
The excitations can be classified according to the num- 
ber of electrons involved in the transitions, leading for 
the 2-electron quantum well to the distinction between 
single- (SE) and double-excitations (DE). In particular, 
the double excitations are of high interest in the dynam- 
ics of correlated electronic systems. Obviously, they can- 
not be captured by mean-field type approaches such as 
Hartree-Fock. For A* = 1, the lowest excitation from the 
ground state, denoted by "a", is a SE of energy ui® D — 
0.61 E*. It is followed by a DE('T'), with wf D = 2.17 E^, 
another SE ("b"), at uj^ d = 2.95 Eq, two DEs, ("2") at 
ujED = 4 54 E * r respectively, ("3") at wf D = 5.31 E* a and 
another SE ("c") with energy D = 6.89 Eq. 

Consider now the results from the GKBA calculations. 
It can be clearly seen from Fig. [2] that it shows single 



excitations as well as double excitations, although the 
quality of their description differs significantly. While 
the SEs are very well described and are practically iden- 
tical to the exact results over the whole range of cou- 
pling parameters A*, the double excitation energies only 
coincide in the limit A* — s- 0, and the deviations from 
the exact result increase approximately linearly with A*, 
leading to a relative error of the order of 20 % at A* = 2.5 
for the lowest DE ("1"). Thus the primary conclusion 
is that our NEGF approach within the GKBA is indeed 
able to reproduce the double excitations in the spectrum. 
At the same time their energy appears with an incorrect 
A* -dependence in the present second Born approxima- 
tion for the self-energy suggesting that not all processes 
leading to these excitations are captured. For complete- 
ness, we note that, for coupling strengths A* > 1, addi- 
tional excitations appear in the spectrum that cannot be 
attributed to real excitations in the system (cf., e.g. "E"). 



To analyze the numerical performance of the GKBA 
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we tested it against full two-time calculations, both us- 
ing the same second order Born self-energy. Fig. [3] 
shows the scaling of the computation time with the 
propagation time T on a single standard CPU. From the 
graph it can be seen, that the full propagation scales with 
T to a power of greater than 3, we expect that it will 
converge to a scaling of oc T 3 for longer time propaga- 
tion. The GKBA, in contrast, scales only as T 2 . This 
figure shows that the GKBA allows to increase the prop- 
agation duration T by three or more orders of magni- 
tude compared to two-time calculations. This, in addi- 
tion to the significant reduction of memory consump- 
tion, paves the way to much longer propagation times 
in the future. 



V. CONCLUSIONS AND OUTLOOK 

In this contribution we have shown recalled the idea 
and previous applications of the generalized Kadanoff- 
Baym ansatz and we demonstrate how to extend it to 
inhomogeneous, finite systems. A key for an efficient 
and consistent simulation was to correctly provide a 
correlated initial state via adiabatic switching of the in- 
teraction. For a first test of the accuracy of the ap- 
proximations we have applied the formalism to a 2- 
electron quantum well model system and have studied 
the ground state dipole excitation spectrum for different 
coupling parameters. While our approach is easily ap- 
plicable to systems containing more particles, the two- 
electron case allows for a benchmark against exact diag- 
onalization results. 

Our numerical results confirmed that the GKBA cor- 
rectly recovers, besides single-particle excitations, also 
double excitations. These are presently of high inter- 
est for many applications in semiconductor optics and 
transport but cannot be obtained by standard tools such 
as time-dependent Hartree-Fock. At the same time, we 
have found that while the SE are reproduced with high 
accuracy, the DE are correctly captured only at small 
coupling. Inherent to the GKBA in second Born approx- 
imation is an incorrect coupling parameter dependence 



of the DE energy of the order oc (A*) , in contrast to 
the correct scaling of oc (A*)~ 2 . This is in good agree- 
ment with our similar findings for the 4-electron quan- 
tum well case |28]. Since a similar scaling is observed 
for the full two-time propagation within the second or- 
der Born self-energy we conclude, that this behavior of 
the DE energies is not a deficiency of the GKBA but in- 
dicates the limitations of the involved (weak coupling) 
second Born approximation for the self-energy. Evi- 
dently, higher terms in the Born series are required to 
restore the correct scaling. Therefore, in future work we 
will study higher order approximations for the selfen- 
ergy such as T-matrix- or GW-approximation. It will 
be interesting to see whether the GKBA performs simi- 
larly well allowing again to omit the complicated inte- 
gral terms in the full equations 1 23 24} . 

Our results based on the GKBA open up a broad vari- 
ety of new many-body applications of inhomogeneous 
finite systems. In our recent work |2j |27| we demon- 
strated that these systems become tractable by using the 
FEDVR representation. Still there were essential limita- 
tions of full two-time calculations in terms of computa- 
tion time and memory requirements. These limitations 
can now be mitigated to a large extend with the help of 
the GKBA by increasing the duration of the time prop- 
agation by more than three orders of magnitude. Not 
only does this allow for a much more accurate com- 
putation of electronic spectra, including double excita- 
tions. At the same time, true nonequilibrium problems 
such as nonlinear excitation and relaxation dynamics or 
pump-probe problems in inhomogeneous systems are 
now within reach of NEGF simulations. 
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